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Abstract - The Newns-Anderson model is ubiquitous in studies of the molecular transport in the 
presence of solvent (outer) reorganization. The present work demonstrates that intramolecular 
reorganization can also be significant for the transport through molecules with floppy degrees of 
freedom, for which the Newns-Anderson model can be extended. The expressions of the model 
parameters deduced from electronic structure calculations for (4, 4')-bipyridine (44BPY) quan- 
titatively differ from those characteristic for outer reorganization due to strong intramolecular 
anharmonicities. These expressions can be utilized as input in transport calculations for 44BPY- 
based molecular junctions of experimental interest [Science 301 (2003) 203, J. Am. Chem. Soc. 
130 (2008) 16045, Nano Lett. 12 (2012) 354]. 



' Introduction. — When an electron subjected to a 
, source-drain bias V is transferred from the negative elec- 
fSI ' trode into (say,) the empty LUMO of a neutral molecular 
, system M*' embedded in a molecular junction, a transient 
'radical anion is created in a first process (M° — 7>M*~), 
, which is subsequently transferred to the positive electrode 
' (M*~ — >M°). If these charge transfer processes are fast 
, (strong electrode- molecule coupling), they cannot be con- 
. . ' sidered as separate events; electron transport proceeds by 
^ , coherent tunneling. The nuclei arc frozen at the opti- 
^> ■ mum molecular geometry Qq. This picture is ubiquitously 
5^ , adopted to describe vacuum molecular junctions within 
. NEGF (nonequilibrium Green's functions) approaches [1]. 
For molecular junctions immersed in electrolytes [2], 
this picture has a counterpart called adiabatic transport 
[SHH], which can be summarized as follows. At a given sol- 
vent configuration Q, the molecular device is traversed by 
a so-called partial electric current i{V;Q) [SHS] which is 
the result of coherent tunneling. The slow solvent dynam- 
ics with respect to the electronic motion legitimates the 
usage of a classical effective coordinate Q. Changes in Q 
model solvent's dipoles/charges that rearrange to stabilize 
the extra anion's charge and induce significant variations 



of the Gibbs adiabatic free energy Q — Q{Q\V) |3H5] at 
a a time scale much shorter than the measurement times. 
Therefore, as schematically depicted in fig. [TJ to compare 
with experiments one has to consider the so-called total 
current I{V), which is obtained by ensemble averaging the 
partial current with a weight function V{G{Q)', V) that de- 
pends on g{Q;V) [cf. eq. ^ below]. 

The Newns-Anderson model [Bl[7] is widely employed in 
these studies [31[51IH]. Since early studies on atoms ab- 
sorbed on a metal surface this model continues to be 
used to describe phenomena of recent interest for molec- 
ular electronics, e. g., transition voltage spectroscopy 
[9l-[T2]. It models the molecule as a single level Eq, which 
mediates the tunneling between the source and the drain. 
This level models the closest molecular orbital (LUMO for 
n-type conduction) to the metal's Fermi energy ep- For a 
molecule immersed in electrolytes, its energy Eq = £o(Q) 
fiuctuates due to solvent reorganization [3l[5l[8|. 

The aim of the present Letter is fourfold: (i) to extend 
the Newns-Anderson model to include intramolecular nu- 
clear reorganization; to show that this extension (ii) is rele- 
vant for molecules of current interest in molecular electron- 
ics but (iii) the underlying functional dependencies differ 
from those describing the solvent reorganization; (iv) to 
deduce these dependencies from electronic structure cal- 
culations. Thus, information will be provided that can be 
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Fig. 1: Schematic representation of the adiabatic molecular 
transport in the Newns- Anderson model [eqs. ^ and ((SJ]. 
The total current represents a Q-ensemble average of LUMO- 
mediated tunneling, whose energy £o(Q) fluctuates within an 
energy range ~ 2A [cf Eq. ([H])]]. See the main text for details. 

subsequently used as input for transport studies based on 
realistic parameters. Transport calculations and compari- 
son with experiments [T51ll5[l55] deserve a special analysis 
and will make the object of a separate publication. 

Specific details. — As a concrete molecule, (4, 4')- 
bipyridinc (44BPY) will be considered in the present 
study. In view of its special structure (cf. fig. [2]), with two 
active nitrogen atoms in para positions, 44BPY is par- 
ticularly suitable for simultaneous binding to two metallic 
electrodes. It has been utilized to demonstrate for the first 
time the possibility of repeated formation of molecular 
junctions |13| . Compounds based on 44BPY, commonly 
known as "viologens" , attracted considerable attention for 
many decades. 44BPY molecules have been incorporated 
in redox active tunneling junctions to demonstrate the 
LUMO electrolyte gating [2]. Several theoretical studies 
devoted to electron transport in 44BPY [T5HT8] consid- 
ered the coherent tunneling at fixed geometry but did not 
examine the impact of intramolecular reorganization. 

The quantum chemical calculations for the present 
study have been done with the Gaussian 09 package [H] 
at density functional theory (DFT) level by using the 
B3LYP functional. Basis sets of double-zcta quality aug- 
mented with diffuse functions (Dunning aug-cc-pVDZ) for 
the light atoms and with relativistic core potential (cc- 
pVDZ-PP from ref. [5D]) for gold have been employed. 

The spinless Newns- Anderson model. The cen- 
tral assumption of the transport approaches based on the 
Newns- Anderson model is that electric conduction is dom- 
inated by a single molecular level. As shown below, this 
should certainly be the case for molecular junctions based 
on 44BPY pHTS]. From A-DFT calculations [H], we de- 
duced a HOMO-LUMO gap A = Ec+Ea~2En = 8.5 eV. 
This quantity, expressed in terms of the energies (E) of 
the cation and anion radicals, and neutral species (sub- 



scripts C, A, and N, respectively) is the counterpart of 
the so-called charge gap used in solid state or mesoscopic 
physics (see ref. [22] and citations therein). Screening ef- 
fects narrow down this gap [23, but it certainly exceeds 
the Kohn-Sham HOMO-LUMO gap {Aks = 4.97 eV 
|16jV which is known to drastically underestimate A. 
If the electrodes' Fermi level were located midway be- 
tween HOMO and LUMO (a situation wherein the Newns- 
Anderson model would inherently fail, since both HOMO 
and LUMO should contribute significantly), the trans- 
mission (Gamow formula) T w exp(-2d^2m/n^A/2) 
through an underestimated energy barrier of a height 
A/2 = Aif5/2 = 2.49eV and a spatial extension d = 
^'^N-Au + '44BPY - 2 X 2.4 -f 7.11 A~ 11.9 A would yield 
a conductance G/Gq = T « 10~^, which is completely at 
odds with the experimental values (G/Go ^ 10"^ — 10~^ 
[Hm]). Here, Gq = 2e^/h ^ 77.48 /^S is the conduc- 
tance quantum. To conclude, the assumption of a domi- 
nant MO appears to be reasonable for 44BPY. Whether 
the LUMO (clectron/n-type conduction) or the HOMO 
(hole /p- type conduction) is the dominant MO cannot be 
determined from transport measurements in two-terminal 
setups alone. This issue can be addressed, e. g., in elec- 
trolyte gating [2] or thermopower studies [T5|. Because 
they revealed an n-type conduction, we will restrict our- 
selves below to the case of a LUMO-mediated conduction. 

Within the most general Newns- Anderson model [6j|7], 
the single MO of energy Eq it contains can be empty 
("-t.i = 0)' single {n^ + ni = 1), or doubly (ni- = = 1) 
occupied, corresponding to the neutral, radical anion, 
and dianion species, respectively. The second-quantized 
Hamiltonian has the expression 

H= ^0 (Q) + Un^n^ + £ph (Q) • (1) 

where n„ = c^Ca denote electron number operators. A 
Hubbard-type interaction accounts for the Coulomb repul- 
sion U between the two spin directions. In an STM-setup, 
the molecule is coupled to two electrodes (substrate s and 
tip t), The average molecule-electrode couplings Tg^t de- 
termine a nonvanishing level width characterized by the 
broadening functions Tg^t oc j [T]. 

The active MO is coupled to classical intramolecular 
(and, if the case, solvent) vibrational modes Q that re- 
organize upon charge transfer. They modulate the MO 
energy Eq — £o (Q) and store an energy £pk. (Q). This 
yields a Q-depeiidence of the total energy E (Q) = (H). 
For convenience, the energy at the neutral optimum Qo 
will be taken as energy zero, Ejy (Qo) = 0. 

An important issue in molecular transport is whether 
the double occupancy of the active MO is significant or 
not. Albeit entirely correct only if a single-particle pic- 
ture (on which the DFT description relies) is applicable, 
the analysis can be done by observing that eq. ((U can be 
described in terms of two single electron states of energies 
£i(Q) = £o (Q) and £2 (Q) = £0 (Q) + C/ (ei < £2, U > 0). 
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The charge transfer efficiency is determined by the energy 
differences £1^2 — £f- Fluctuations in £0 (Q) induced by 
phonons are of the order of the reorganization energy A, 
which typically amounts a few tenths of electron volts [5] . 
So, a doubly occupied LUMO gives a significant contri- 
bution only if U is at most of the order of A; otherwise, 
the state of energy £2, too high above the Fermi level, is 
blocked, and only that of energy £1 is relevant. 

The Coulomb blockade parameter U cannot be di- 
rectly determined from transport data in a simple manner 
[24H26]. Within DFT, U ^ Uks obtained from energy 
splitting of the Kohn-Sham LUMO "orbitals" is Uks = 
1.54, 1.52, and 1.26 eV for a 44BPY molecule in vacuum, 
in (aqueous) solution, and in solution with one gold atom 
attached at each of the two nitrogen atoms, respectively. 
Similar to the case of the DFT HOMO-LUMO gap, Uks 
drastically underestimates U . A substantially higher value 
is obtained via the more adequate method of energy dif- 
ferences based on eq. ^,U = En + Ed - 2Ea = 1.92 eV 
(instead of 1.26 eV), for the last of the three aforemen- 
tioned situations (subscript D stands for dianion). Still, 
what is really important for the present purpose is that 
U is much larger than the reorganization energies (see be- 
low). Transfer processes with a doubly occupied LUMO 
are energetically too costly and can be ignored in electron 
transport through 44BPY. So, one can safely employ a 
spinless Newns- Anderson model Hamiltonian (n = c^c) 

H^eo (Q) n + Eph (Q) , (2) 

as done in existing studies, e. g., refs. [3l[5l[8]. 

The total energies of the radical anion Ea{Q) and neu- 
tral species i?Ar(Q) can be used to microscopically com- 
pute the Q-dependence of the parameters entering eq. ([2]) 

£ph[Cl) = Sjv(Q); £o(Q) = Sa(Q) - EN{q,). (3) 

Notice that the LUMO energy £o(Q) expressed above is 
measured with respect to the vacuum. The reorganization 
energies of the radical anion (A^) and the neutral (A^v) 
species are important quantities defined by 

\n = En (Qa) ' En (Qo) ■,Xa = Ea (Qo) - Ea (Qa) , 

(4) 

where Qa denoted the radical anion optimum geometry. 

All the phenomenological approaches of molecular 
transport in electrolytes based on the spinless Newns- 
Anderson model of which we are aware (e. g., Refs. |^ 
[4,5,8;27J) assume harmonic Q-depcndencics of En and 
Ea, and this yields (Q = {Q,^}; henceforth Qo = 0) 

£ph{Q.) = \Y.^uQI, (5) 
feo(Q) = £o(Q)-£o(Q = 0) = -^.g.Q.. (6) 

V 

In the cases where eqs. (O and (|6]) apply, the above reor- 
ganization energies are equal 

A. = Aa = A.e£- 

V 




Fig. 2: The neutral 44BPY molecule consists of two pyridine 
rings twisted by S = 37.2°. 

By linearly joining the minima of the radical anion and 
neutral species, one can introduce an effective coordinate 
7?., Qv = TZgy/uji, and recast Eqs. (O and ([6]) as |27] 

£ph{n) ^XTZ"; eoin) = £0 - 2X11. (8) 

The configurations Qo and Qa correspond to 7^ = and 
TZ = 1, respectively. 

Extending the Newns- Anderson model to de- 
scribe intramolecular relaxation. Most in- 
tramolecular modes are fast (frequencies lJi, comparable 
to or higher than t) and should be treated quantum me- 
chanically [28l[29| . Fingerprints of the associated inelastic 
tunneling are, e. g., the well-known peaks in the second 
derivative (Pl/dV^ at resonant voltage values {eV = huju). 
They negligibly reorganize and therefore are not interest- 
ing in the present context. But if (like in 44BPY, see 
below) slow intramolecular vibrations exist that are suf- 
ficiently strongly coupled to the LUMO, their reorgani- 
zation during electron transfer is significant and should 
be considered. The problem is to investigate whether a 
certain molecule possesses such modes and to scrutinize 
whether they can be described by eqs. ([5|) and ©• 

To this aim, we have performed electronic structure cal- 
culations for 44BPY. The values of the reorganization en- 
ergies found from eqs. (j4|) were found different: \n = 
0.353eV and Aa = 0.224cV. The inequality Ajv ^ Aa 
demonstrates that the inner reorganization of the 44BPY 
molecule cannot be correctly described by eqs. (O and ([5]). 

The reorganization energies Aat^a computed via eqs. (U) 
represent "global" quantities emerging from electronic 
structure calculations, wherein all intramolecular vibra- 
tional modes are included within a classical picture. To 
gain further insight, we have split the reorganization en- 
ergy into contributions of individual molecular vibrational 
modes. Composed of 20 atoms, the nonlinear 44BPY 
molecule has 54 normal vibrations. In the D2 point group 
symmetry, they are distributed as 14^1 4- 12Bi + 14i?2 + 
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14i?3; lOA + 9Bi + 9i?2 + 8-63 are in-plane vibrations and 
AA + 3Bi + 5B2 + 5i?3 are out-of-plane vibrations [3D]. 
Our extensive calculations confirmed the expectation that 
significant contributions to the reorganization energy arise 
from the in-plane normal modes with A symmetry. Out of 
the ten in-plane modes of A symmetry we have identified 
two modes that dominate the inner reorganization. As ex- 
pected, they are related to the main structural differences 
between 44BPY0 and 44BPY*-. 

One of these modes (normal coordinate Q4e) is re- 
lated to the so-called quinoidal distortion [30] of the neu- 
tral molecule upon electron attachment (44BPY° 4- e~ — > 
44BPY*^), with a shortening of the inter-ring C-C bond 
and of the C-C bond parallel to it, and a lengthening 
of the C-C bond between them as well as of the C-N 
bond. Adiabatic energy curves EN^AiQio) reveal a virtu- 
ally perfect harmonic behavior, which agree with eqs. ([5]) 
and (|6]) and yield equal partial reorganization energies 
A^^'' = A^^''. Because its frequency (the computed value 
LOiQ — 1642 cm" ^ ~ 0.2 eV excellently agrees with the 
strong Raman band observed experimentally at 1645 cm~^ 
[50] ) exceeds typical r^^t-values [TUHT^ . this mode is too 
fast to significantly reorganize by the classical thermal ac- 
tivation considered in this study. 

The inequality \n ^ traces back to the most ap- 
pealing feature of the molecular structure, namely the 
inter-ring twisting angle 6. The mode directly related 
to the inter-ring torsional motion represents the floppy 
(label /) degree of freedom of 44BPY. Adiabatic energy 
curves i?jv,A(Q/) along the normal coordinate Qf exhibit 
strong anharmonicities. Its (harmonic) frequency is very 
low [ujf ~ 62cm"^ in the isolated molecule, cf. table 
d]), but the impact on reorganization is important be- 
cause large amplitude oscillations are significant. The 

i f) 

partial reorganization energy of the radical anion \^ ~ 

0.16 eV is almost two times larger than that of the neu- 

( f) 

tral molecule (A)^ ~ 0.09 eV). From the adiabatic En^a- 
curves obtained by quantum chemical calculations, the 
functional dependencies of the model parameters Eph (Qf) 
and EoiQf) can be computed from eqs. ([3]). The results 
of these calculations are collected in fig. [3K- They show 
significant deviations from the linear and quadratic Qf- 
dependencies of £0 (Q/) and £ph{Qf), respectively. The 
latter approximations are represented as dashed lines in 
fig. [3^. Fig. |3] also shows that fourth-order polynomials 

£a{Qf) = -EAy+Uf{eiQf + e2Q} + e3Qj + e4Qj), 
SphiQf) = u:f{Q}/2 + hQ} + UQj). (9) 

very accurately fit the DFT-curves. Values of the vertical 
electron affinity EAy = — £o(0) = i?jv(Qo) - ^a(Qo) and 
the other quantities entering eq. (jH]) are given in table [T] 

Beyond the case of an isolated molecule. 

The results discussed above refer to an isolated 44BPY 
molecule in vacuum. Cases of breakdown of the harmonic 
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Fig. 3: The model parameters eo and £ph of a 44BPY molecule 
(a) in vacuum, (b) in aqueous solution, and (c) with two gold 
atoms attached. The relevant range of the normal coordinate 
of the floppy degree of freedom Q / includes values around the 
minima at Q'j and Qf of the neutral and radical anion species, 
respectively. The results of the DFT calculations (solid lines) 
can be excellently fitted with polynomials [eqs. (|9]) and table 
[T] (represented by points), but they significantly depart from 
the linear and quadratic approximations of eqs. (jSJl and ([5]), 
respectively (dashed lines). See the main text for details. 

approximation for low frequency vibrations in isolated 
molecules are well known in molecular physics [31] . How- 
ever, we are not aware of studies pointing out the failure 
of the harmonic approximation for floppy molecules used 
to fabricate molecular junctions, like in the case discussed 
above. For molecular transport, it is also important to 
consider the case of a 44BPY molecule immersed in an 
electrolyte and linked to metallic electrodes j2l [T3l[T5] . 
In the present study on 44BPY in aqueous solution (la- 
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/3 


h 


ei 


62 


63 


64 


EA, 


44BPY in vacuum 61.9 


0.0399 


0.10038 


-9.0236 


-0.9194 


0.04239 


-0.0007 


0.444 


44BPY in solution 58.18 


0.0634 


0.1144 


-9.6301 


-1.0751 


0.04438 


0.00115 


2.206 


Au-44BPY-Au in vacuum 63.2 


-0.0040 


0.0943 


-9.0110 


-0.2416 


-0.3563 


0.05714 


1.586 



Table 1: Results for the 44BPY molecule in vacuum, in aqueous solution, and with two gold atoms connected to the two 
nitrogen atoms. The units for the frequency of the floppy mode ujf and the vertical electron affinities EAv are cm~^ and eV, 
respectively. The dimensionless coefficients /3,4 and 61,2,3,4 enter the interpolation formulas of eqs. 



bel sol), the solvent has been considered within the po- 
larized continuum model using the integral equation for- 
malism (keyword SCRF^IEFPCM in Gaussian 09). The 
results (fig.|3)3 and table [T|) indicate that the most impor- 
tant solvent effect is an almost constant shift in the LUMO 
energy. It is related to the strong interaction of the an- 
ion with the water molecules. The largest contribution to 
the vertical electron affinity in solution comes from the an- 
ion's solvation free energy AGa = E'^^{CI''/)-Ea{CIa) ^ 
— 2.01eV. Here, Q^°' denotes the optimized geometry of 
the radical anion in solution. 

As a preliminary step in investigating the impact of 
electrodes, we have also considered the case of a 44BPY 
molecule with two gold atoms attached to the nitrogen 
atoms. The corresponding results (fig. [3};) reveal an effect 
qualitatively similar to that of the solvent. 

So, apart from a nearly constant shift of the LUMO 
energy, for Q /-values of interest (cf. fig. [3]) , the model pa- 
rameters EoiQf) and SphiQf) vary within ranges, which 
appear to be little sensitive to the presence of solvents 
or electrodes. The physics behind the similarity exhib- 
ited by the these numerical results is the following. The 
twisting angle 6 of the 44BPY molecule is determined by 
the competition between the 7r-clcctronic interaction of the 
pyridyl fragments, which favors 7r-electrons delocalizcd be- 
tween coplanar pyridyl rings (0 = 0), and the steric repul- 
sion between the ortho H atoms, which is diminished by 
a twisted conformation {6 ^ 0) [32]. The latter prevails 
in the neutral species 44BPY'' (empty /oxidized LUMO), 
which is nonplanar (fig. [2|). By adding an extra electron 
(occupied/reduced LUMO), the energy gain resulting from 
TT-electron delocalization between the two rings outweighs 
the steric repulsion, and the radical anion (44BPY*~) be- 
comes planar. Therefore, a structural transition from the 
twisted to the planar conformation can only be suppressed 
if a significant negative charge is transferred to the 44BPY 
unit. This cannot be achieved by immersing in solution, 
nor even by attaching gold atoms. In the latter case, the 
gold atoms acquire a small negative charge — 0.096e), 
which has a negligible impact on the torsion angle 9. 

Summary and outlook. — The present results em- 
phasize the need to consider the inner relaxation in junc- 
tions based on molecules with floppy degrees of freedom; 
the reorganization energies deduced here (A ~ 0.1 — 0.2 eV) 
are comparable with those for outer (solvent) reorgani- 
zation [2]. To avoid misunderstandings three comments 
are in order, (i) The present analysis of the inner reor- 



ganization (conformational distortions) has implicitly as- 
sumed a given (e. g., atop, hollow, bridge) contact ge- 
ometry; variations in the binding geometry may be larger 
(e. g., refs. but they are of a different nature, (ii) 

The Q f dependence of the partial current and conduc- 
tance discussed here refers to a given molecule (44BPY); 
this is qualitatively different, e. g., from the scaling G oc 
cos^ 8 discussed previously |33fl35j , which refers to various 
derivatives of a molecular family (e. g., biphenyls) wherein 
the torsion angle 9 is varied, e. g., by a bridging (alkyl) 
chain, (iii) 9 cannot be specified in terms of the single 
normal coordinate Q j only. Thermal averaging over 9 (as 
done for biphenils [33]) amounts to consider that, concomi- 
tant with the floppy mode w/, other modes of much higher 
frequency are also thermally activated. Because such fluc- 
tuations are energetically costly in biphenyls, their effect 
is reduced [33j . 

While revealing that a spinless Newns- Anderson is jus- 
tified for 44BPY, the present study has shown that the 
Q-dependencies of eqs. ([5]) and © assumed to work in 
electrolytes are inappropriate for the reorganization of low 
frequency intramolecular vibrations, which are character- 
ized by a pronounced anharmonic behavior. 

New formulas for eo{Qf) and 8ph{Qj) have been de- 
duced from microscopic calculations [cf. eqs. ©], which 
replace the expressions of eqs. ([S]) and (jS]) used for solvent 
reorganization. The Q ^-dependence of the model parame- 
ters is important because the ensemble average needed to 
compute the experimentally relevant quantity, namely the 
total current I{V), requires an integration over Qj 

I{V)^{]{V-Qf))^ I j{V;Qf)P{Qf;V)dQf. (10) 

The fact that 44BPY possesses a single low frequency vi- 
brational mode that significantly reorganizes represents an 
enormous simplification; otherwise, an ensemble average 
involving a Q-integration over all of the 54 internal degrees 
of freedom would be a formidable numerical challenge. 

The thermal weight function V{Qf]V) oc 
exp [— fj((5/; T^)/(/cbT)] is determined by the Gibbs 
free energy G{Qf]V) of the (partial) oxidation state of 
the molecule (i. e., LUMO occupancy < n{Qf]V) < 1) 
linked to biased electrodes. Q depends on Q/, V, and 
(if applicable) on electrolyte's overpotcntial [2j]5]. The 
need to carry out an ensemble averaging is an essential 
aspect, which renders purely ab initio approaches (as used 
for rigid molecules in vacuum) inapplicable; therefore, 
feasible approaches to date have to resort to models. 
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The above results can (and will) be used in subsequent 
studies to deduce the current within the adiabatic trans- 
port approach, for which eq. (jlOp can serve as starting 
point. Expressions for the partial current j{V] Sq^Q f)) 
and the level (LUMO) occupancy n{eo{Qf); V) are known 
[SJISlinnilSZ]- Along with the expression of SphiQf), the 
LUMO occupancy is required to express the Qf- (and V- 
)dependent Gibbs adiabatic free energy Q. One should em- 
phasize at this point that G{Qf', V) is not simply related 
to the adiabatic energies Em,a{Q f)] the former quantity 
characterizes a molecule linked to biased electrodes, while 
the latter quantities pertain to an isolated molecule. 

Earlier studies demonstrated that even a single reorga- 
nizable harmonic mode, eq. ([5]), yields highly nontrivial 
adiabatic ^7-surfaces [5l[8]. In view of the significant an- 
harmonicities embodied in the expressions adequate for 
a floppy degree of freedom, one can expect more com- 
plex adiabatic surface topologies, with the need to distin- 
guish between several important limiting cases. To sim- 
ply motivate this, one should note that, unlike in elec- 
trolytes, reorganization effects in floppy molecules cannot 
be merely characterized with the aid of a single quantity 
(A = Aat = \a)- 

As shown recently [ITJ[T2|, the impact of fluctuations 
in the MO-energy offset eo = eq — ep on the ohmic con- 
ductance can be significant. Molecule-electrode interac- 
tions can be important sources of such fluctuations 
In view of the present study, in molecular junctions like 
those based on 44BPY fabricated experimentally Pll3ll5j , 
one can expect that, even in the absence of other effects, 
reorganization effects represent an important source for 
large eo-fluctuations (Jeq 0.2 — 0.4 eV). As a straight- 
forward application, one can employ eq. (|10p to compute 
conductance histograms, which can be directly compared 
with those available from experimental studies [13[ll5j. 

* * * 
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